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Using the spin-rotation-invariant Green's function method we calculate the thermodynamic quan- 
tities (correlation functions (SoSr), uniform static spin susceptibility x> correlation length £, and 
specific heat CV) of the two-dimensional spin-1/2 J\-J% Heisenberg ferromagnet for J2 < Jf ~ 
0.44|Ji|, where J£ is the critical frustrating antiferromagnetic next-nearest neighbor coupling at 
which the ferromagnetic ground state gives way for a ground-state phase with zero magnetization. 
Examining the low-temperature behavior of \ an d £, in the limit T — > both quantities diverge 
exponentially, i.e., \ oc exp(fe/T) and £ oc exp(6/2T), respectively. We find a linear decrease in 
the coefficient b with increasing frustration according to b — — 5 (Ji + 2J2), i.e., the exponential 
divergence of x an d £ i s present up to Jf. Furthermore, we find an additional low-temperature 
maximum in the specific heat when approaching the critical point, J2 — ► J%. 

PACS numbers: 



I. INTRODUCTION 

Low-dimensional quantum magnets have attracted 
much attention during the last decades*^ They are pre- 
destined to study the influence of strong thermal and 
quantum fluctuations. Much attention has been paid 
to the theoretical investigation of the two-dimensional 
(2D) spin-1/2 J1-J2 quantum Heisenberg antiferromag- 
nct which may serve as a canonical model to study 
the interplay of frustration effects and quantum fluc- 
tuations (see, e.g., Ref. 0). An additional motivation 
to study this model comes from the experimental side^ 
Very recently, several quasi-2D magnetic materials with 
a ferromagnetic nearest-neighbor (NN) coupling J 4 < 
and a frustrating antiferromagnetic next-nearest neigh- 
bor (NNN) coupling J 2 > have been investigated ex- 
perimentally, e.g., Pb 2 VO(P0 4 ) 2 r- (CuCl)LaNb 2 7 ^ 
SrZnVO(P0 4 ) 2 ^^ and BaCdVO(P0 4 ) 2 £i&±2 The 
quite large frustrating J 2 drives these materials out of 
the ferromagnetic phase. The experimental findings have 
stimulated several theoretical studies of the ground-state 
and thermodynamic properties of the J 4 - J 2 model with 
Ji < and frustrating J 2 > 0J£~— It was found that 
the ferromagnetic ground state for the spin-1/2 model 
breaks down at J 2 = Jf « 0.4| Ji| Note that 

for the classical model (spin s —> 00) the corresponding 
transition point is at J 2 = 0.5 1 J±\. 

On the other hand, some materials considered as 
2D spin-1/2 ferromagnets, such as K 2 CuF 4 , Cs 2 CuF 4 , 
Cs 2 AgF 4 , La 2 BaCu0 5 , and Rb 2 CrCl 4 ^— might have 
also a weak frustrating NNN interaction J 2 < Jf . In 
contrast to the previous investigations of the 2D J\- J 2 
model with fcrrromagnetic J^ -18 : 21 which have consid- 
ered predominantly the case of strong frustration, in the 



present paper we will focus on the region of weak frus- 
tration J 2 < Jf . 

Although for J 2 < Jf the ground state remains fer- 
romagnetic, the frustrating J 2 may influence the ther- 
modynamics substantially, in particular near a zero- 
temperature transition. This has been demonstrated 
for the one-dimensional (ID) frustrated ferromagnet in 
Ref. where a change in the low-temperature behav- 
ior of the susceptibility and the correlation length as well 
as an additional low-temperature maximum in the spe- 
cific heat have been found when approaching the zero- 
temperature critical point from the ferromagnetic side. 

The Hamiltonian of the system considered in this paper 
is given by 



H = Ji s * s j + J2 s * s ^ ; ^1 < ; J 2 > 



(1) 



where (S^) 2 = 3/4, and denotes NN and 

denotes NNN bonds. The unfrustrated 2D ferromag- 
net (J 2 = 0) has been widely investigated, e.g., by 
the modified spin-wave theory^ renormalization group 
approaches) 30 ' 31 the quantum Monte Carlo method^ - — 
and by a spin-rotation-invariant second-order Green's 
function method (RGM)^— However, in presence of 
frustration (J 2 > 0) the choice of appropriate methods 
for studying temperature dependent quantities of this 
system is restricted due to frustration and dimension- 
ality. For instance, the powerful quantum Monte Carlo 
method is not applicable due to the minus sign problem, 
whereas finite-temperature density matrix renormaliza- 
tion group studies are restricted to ID systems. 

So far, only a few theoretical papers deal with the ther- 
modynamics of the model for J 2 > 0, see, e.g., Refs.dll 
and where finite lattices of N = 16 and N = 20 
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sites arc considered. As we will sec below, for these lat- 
tice sizes the finite-size effects at low temperatures are 
large. Therefore, methods studying infinite systems are 
highly desirable. Among others, the Green's function 
method is a powerful tool to study mag netic systems 
at finite temperatures, see e.g. Refs. l40M42l and refer- 
ences therein. A specific variant appropriate for frus- 
tated systems is the RGM. In the present paper, we use 
the RGM to investigate thermodynamic properties of the 
frustrated model ([T]), focusing on J 2 < Jf ~ 0.4|Ji|. The 
RGM has been applied successfully to low-dimensional 
(frustrated) infinite quantum spin systemsi 28 ' 35-39 ' 43 " — 
In particular, for the ID spin-1/2 Heisenberg ferromag- 
nct it was shown that the RGM reproduces Bethe-ansatz 
results . 28 ' 37 Since the RGM can be formulated also for 
finite systems, we use full exact diagonalization (ED) of 
finite square lattices of N = 16 and N = 20 sites to 
compare the RGM with ED data for finite lattices. 

The paper is organized as follows: In Sec. QI] the RGM 
applied to model ([1} is presented. The results for the 
thermodynamic quantities are discussed in Sec. IIIII Fi- 
nally, a summary is given in Sec. IIVI 
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where A M = 1 + 2a (2C M + C 2 , - 5Ci, ), K 2 , 2 = 1 + 
2a (2C 2 ,o + C 2 ,2 - 5Ci,i), A\ 2 = 4a (<7 1)2 - C 1>0 ), and 
A'2,1 = 4a (Ci j0 + Ci, 2 - 2C14). From the Green's func- 
tion ([|5]) the correlation functions Cr = ^ q C q e lqR 
of arbitrary range R are determined by the spectral 
theorem 
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(8) 



where n(oj q ) = (e^< l / T — l) 1 is the Bose function. Tak- 
ing the on-site correlator Cr = o and using the oper- 
ator identity Sf S~ = \ + Sf we get the sum rule 
^0 = ^XiqCq = 5" ^ ne un iform static spin suscep- 
tibility x = hniq^oXq: where x q = Xq(^ = 0) and 
Xq (<*>) = \xt~ M, is given by 



II. ROTATION-INVARIANT GREEN'S 
FUNCTION METHOD 



A 



J2 kJ k Ci,k-i, A = Yj kJ kJiKkj- 



(9) 
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To calculate the thermodynamic quantities of the 
model (J]), we use the RGM^L 3 -^^ for determin- 
ing the transverse dynamic spin susceptibility Xq ~ ( w ) = 
— ((5+ ; SZq))u, where ((...;.. denotes the two-time 
commutator Green's function^ Taking the equations of 
motion up to the second step and supposing spin rota- 
tional symmetry, i.e., (Sf ) = 0, we obtain 



w 2 ((S+;Sl q >>. =M,+ <<-S+;Sl q )). , 

M q = ([[5+, J ff],5_ q ]) , -§+=[[S+,H],H] 
The moment M q is given by the exact expression 
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E JkCi,k-i (1 



fc=l,2 



(2) 



(3) 



(4) 



where C n , m = C R = (S+S R ) = 2(5 Z 5 R ), R = ne x + 
me y , 7 q 1 - ) = (cosffc +cosq y ) /2, and 7 q 2 ' 1 = cos q x cos q v . 
The second derivative — 5+ is approximated in the spirit 
of the schemes employed in Refs. mHHHiifi! That 
is, in —Sf we adopt the decoupling 

SfSfSZ = a itk {SfSz)Sf + a jtk (SfSz)Sf. (5) 



Following the arguments of Rcf. |28j we put a,^ = a in 



the whole temperature region. We obtain — 5+ = u>„S^ 



and 



X q (w) = -«5+;5_ q )) (1 



(6) 



The correlation length may be calculated by expanding 



Xq around q 
find 



0, 



28.35,36 



Xq = x/ [1 + e 2 {ql+ql)\- We 



2a(J 1 2 C 1 , + 2Ji J 2 (Ci, + Ci,i) + 4jf Ci,i) 



(10) 



The ferromagnetic long-range order, occurring at T = 
only, is reflected by a non- vanishing quantity C (see, e.. 



Ref. [36|) according to C R (0) = ^ £ q#0 C q (0)e i i R + C, 



where C, typically called condensation term^ is con- 
nected with the magnetization m by m 2 = 3C/2. Due 



_ M q (0) 



must 



to the exact result Cr^o(0) = g, C q (0) — 9^~(o) 
be independent of q. Examining Eqs. fU) and ([7]), this 
requires Afc^O) = which leads to a(0) = | and 
C q (0) = |. Hence, in case of a ferromagnetic ground 
state we find at zero temperature w q = 2p s q 2 (|q| << 1), 
~~ j(Ji + 2^2) (where p s is the spin-stiffness) and 



Note that the sum rule Cq = h is 



Ps = -4W1 

Cr(o) = i<y F 

fulfilled. In Eqs. © and JTDJ}, we have A(0) = so that 
X and £ diverge as T — >• indicating the ferromagnetic 
phase transition. 

To estimate the zero-temperature transition point J 2 
(see Sec. IIII A[) . we determine the ground-state corre- 
lation functions and the ground-state energy per site, 
E = i{J\C\fl + J 2 Ci,i), for J 2 > J|. Thereby, we take 
into consideration the possible existence of a quantum 
collinear phase described by the magnetic ordering vector 
Q = (n, 0) or (0, it) and proceed along the lines indicated 
in Ref. Hj. 

Since we will compare the RGM data with ED data for 
finite lattices, see Sec. IIIII we have to adopt the RGM to 



3 



finite N. In this case the quantity C is not related to the 
magnetization and stays nonzero in the whole tempera- 
ture region. In Ref. |45| it was shown that C = 2Tx/N. 

To evaluate the thermodynamic quantities, a coupled 
system of six (for finite systems: seven) nonlinear alge- 
braic self-consistency equations, including the sum rule 
Co = 5, has to be solved numerically to determine the 
correlators C^o, Ci,i, C^o, C2,i, C2,2> (C) and the ver- 
tex parameter a. To solve this system we use Broyden's 
method^ which yields the solutions with a relative error 
of about 10 -8 on the average. The momentum integrals 
are done by Gaussian integration. To find the numerical 
solution of the RGM equations for T > 0, we start at 
high temperatures and decrease T in small steps. Below 
a certain (low) temperature Tq(J2) no solutions of the 
RGM equations (except at T = 0) could be found, since 
the quantity A(T, J2) in Eqs. (J5J and (|10l) becomes ex- 
ponentially small which leads to numerical instabilities. 

Presenting our results in the next section we put 3\ = 
-1. 



find solutions of the RGM ground-state equations down 
to J 2 = 0.46. This value lies certainly above the transi- 
tion point Jf . To fix Jf we may extrapolate the ground- 
state energy obtained for case (ii). The crossing point 
of both energies at J2 = 0.44 can be considered as the 
RGM estimate of Jf , see the inset of Fig. [TJ This value 
is in reasonable agreement with values for Jf reported in 
other P apersi±2iI2iI2i2ii22 

The behavior of the spin-spin correlation functions at 
T = near the transition point Jf i s illustrated in Fig. [T] 
For J2 > Jf there is a noticeable variation in the cor- 
relation functions. With increasing J2, at J2 > 1, the 
spin-spin correlation functions approach the correspond- 
ing values of the limiting case J2 3> 1. This behavior 
is in qualitative agreement with data from Lanczos ED 
and coupled-cluster method^ At the critical point the 
correlation functions jump to the exact results of the fer- 
romagnetic phase, i.e., (SoSr) = 0.25. Together with 
the kink in the energy, this indicates a first-order tran- 
sition. These results confirm the findings in previous 
papers.— iI2i22 



III. RESULTS 



Phase transition at T = 



In a first step we use the RGM as described above to 
determine the transition point Jf , where the ferromag- 
netic ground state gives way for a ground-state phase 
with zero magnetization. For that we solve the RGM 
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B. Thermodynamic properties 

Now we investigate the thermodynamic properties of 
the model for J2 < Jf ~ 0.44, where the ground state 
is ferromagnetic. First we consider the temperature de- 
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FIG. 1: Ground-state spin-spin correlation functions (SoSr) 
in dependence on frustration J2. (Note that for J2 < J| ~ 
0.44 all (SoSr) coincide.) Inset: Ground-state energy E ver- 
sus J2 (solid line: RGM data, dashed line: extrapolated RGM 
data to fix the transition point, see also text). 



equations at T = (i) starting from J2 = with increas- 
ing J 2 and (ii) starting from large J 2 S> 1 with decreasing 
J2. In case (i) the RGM equations at T = can be solved 
until the classical transition at J2 = 0.5. In case (ii) we 



FIG. 2: NN and NNN spin-spin correlation functions for J2 = 
0, 0.15, and 0.3, from top to bottom, calculated by RGM for 
N —>■ 00 (lines) and ED for N = 20 (filled symbols). 

pendence of the spin-spin correlation functions. The NN 
and NNN correlators are shown in Fig. [5] for various 
values of J2. For comparison we show also ED data for 
N = 20 sites. The RGM results agree qualitatively with 
the ED data. With increasing frustration, the correla- 
tion functions decrease more rapidly with temperature. 
Interestingly, for large frustration the NNN correlation 
function changes its sign at a certain temperature, e.g., 
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for J 2 = 0.3, at T w 1.25 (RGM data). The faster decay 
of the spin-spin correlators due to frustration is related 
to the decrease in the correlation length £ with increasing 
J 2 , see the discussion below. 



temperature both \ and £ decrease rapidly with J2 
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FIG. 3: Uniform static spin susceptibility x for J2 = 
0,0.1,0.2,0.3,0.4, from right to left. The inset shows the 
susceptibility for T = 0.3 in dependence on the frustration 
J2. The solid line represents an exponential fit, x(T = 0.3) fs 
132.8 exp(— 9.7J2), of the data points (filled squares). 




FIG. 4: Correlation length £ for J 2 = 0,0.1,0.2,0.3,0.4, 
from right to left. The inset shows the correlation length for 
T = 0.3 in dependence on the frustration J2. The solid line 
represents an exponential fit, £(T = 0.3) w 13.9exp(— 6.U2), 
of the data points (filled squares). 

The uniform static spin susceptibility \ an d the cor- 
relation length £ depicted in Figs. [3] and |4] show a sim- 
ilar behavior. Due to the ferromagnetic ground state, 
both quantities diverge at T = exponentially, see be- 
low. With increasing frustration J 2 the rapid increase 
in both quantities is shifted to lower temperatures. As 
shown in the insets of Figs. [3] and SI at a certain fixed 
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FIG. 5: Uniform static spin susceptibility x f° r ( a ) J2 = 0.1 
and (b) J 2 = 0.3 for N = 16, 20, 64, 400, thermodynamic limit 
(lines, RGM) and N = 16,20 (filled symbols, ED), from left 
to right. Note that for J2 = 0.3 the curves for N = 400 and 
TV — !> 00 almost coincide. 

Since the susceptibility (together with the correlation 
length) is an important quantity to analyze the critical 
properties for T — > 0, we first test the quality of our 
RGM approach in more detail by comparing the RGM 
data for N = 16 and N = 20 with ED data for the same 
system sizes, see Fig. \5\ It is obvious, that the RGM 
and the ED data for x are in excellent agreement. This 
is consistent with our previous finding^ that for the ID 
Heisenberg ferromagnct the RGM data for x(T — > 0) 
and £(T — > 0) are in perfect agreement with exact Bethe 
ansatz results. Moreover, the RGM calculations of x for 
finite N allow to estimate the magnitude of the finite-size 
effects. As can be clearly seen in Fig. [51 the finite-size 
effects at T < 0.6 (T < 0.4) for J 2 = 0.1 (J 2 = 0.3) be- 
come very large. Since the correlation length £ becomes 
smaller with increasing J 2 (see the inset in Fig. [4]) , the 
finite-size effects become less pronounced with growing 
frustration. Nevertheless, based on our data we argue 
that the finite systems of N = 16 and N = 20, accessible 



5 



by exact full diagonalization, are not representative for 
the thermodynamic limit, see also Refs. [l6| and[l3- Note 
that for the ID case, the ED data for N — 20 are in 
good agreement with RGM data for N — > oo»2& Consis- 
tent to that observation, in two dimensions the system of 
N = 400 sites with the same linear extension is already 
close to the thermodynamic limit, see Fig. [5] 

Next we use the RGM for TV — > oo to investigate 
the critical behavior of x an d £ for T — > in more 
detail. We start with a brief discussion of the low- 
temperature behavior of the unfrustrated fcrromagnct. 
Contrary to the ID case, where for J 2 = exact Bcthe 
results arc available j^ 7 - we have only approximate re- 
sults for the 2D model. From low-temperature expan- 
sions of the susceptibility and the correlation length for 
the model with J2 = using rcnormalization group 
approaches the modified spin-wave threory^ and 
the RGM 3 - 7 . it is known that for T -> the susceptibility 
behaves as x ^ T s exp (b/T) and the correlation length 
as £ oc T CT exp ((3/T). While different leading exponents 
s and a for the (less important) preexponential factor 
were obtained by different methods, the exponential di- 
vergence is obtained by all authors^ -31 ' 37 However, dif- 
ferent values for the coefficients 6 and fi were reported, 
namely, b(J 2 = 0) = tt/2 using RGM^l b{ J 2 = 0) = tt us- 
ing modified spin- wave theory^ 9 - or renormalization group 
approach^ and b(J 2 = 0) = 0.132 77T using a different 
version of the renormalization-group approach^ In all 
these papers ^ 31 ' 37 it was found that the coefficients in 
the exponents fulfill the relation j3 = 6/2. We mention 
further, that early numerical studies based on quantum 
Monte-Carlo calculations (using, however, x an d £ data 
only for quite large temperatures) give 6( J 2 = 0)j~ 4.5 « 
1.43tt (Ref. H| and 0(J 2 = 0) = 0.254tt (Ref-H). Let 
us finally argue, that the results for the coefficient b ob- 
tained by modified spin-wave theory^ 9 - as well as renor- 
malization group approach 3 -^ and by the RGM 3 - 7 - seem 
to be most reliable, since these methods are well-tested. 
Moreover, for the ID spin- 1/2 Hcisenberg ferromagnct 
it was shown that the RG M 28 i 37 as well as the modified 
spin- wave theory^ 9 - reproduce the exact Bethe-ansatz re- 
sults for the low-temperature behavior of x an d £■ Hence, 
it is to some extent surprising, that there is a difference 
by a factor of 2 in the coefficient b between the value ob- 
tained by modified spin-wave theory^ 9 - (as well as renor- 
malization group approach 3 -^) and the RGM for the 2D 
unfrustrated ferromagnet. This discrepancy is a known 
but unresolved problemi 35 i 37 However, we believe that 
our general conclusions concerning the exponential diver- 
gence of the frustrated model, see the discussion below, 
are not affected by this problem. 

Next we use the RGM results to determine the co- 
efficients b and j3 as functions of J 2 . We assume that 
the general low-T behavior of these quantities is pre- 
served for < J2 < Jf, cf. Ref. HH, and fit our nu- 
merical RGM data for x an d £ &t low temperatures us- 
ing the ansatzes x = (o-oT _1 + a\ + a 2 T) cxp (y) and 



£ = ^(aoT- 1 + ai + a 2 T) exp UjA. Note that the lead- 
ing power T -1 in the preexponential functions was de- 
rived for the unfrustated case with the RGM in Ref. l37l . 
We consider 6, ao, 01, and a 2 as well as /3, ctQ, a±, and a 2 
as independent fit parameters. Recall that we can calcu- 
late RGM data only for temperatures down to a certain 
T , where T (e.g., T = 0.161 for J 2 = and T = 0.052 
for J 2 = 0.4) is the lowest temperature where the sys- 
tem of RGM equations converges (see Sec. [TTJ) . Thus, 
in addition to the leading power T^ 1 in the preexpo- 
nential function it is reasonable to consider higher order 
terms to achieve an optimal fit of the RGM data. For 
the fit we use 500 equidistant data points in the interval 
T ... T + T cut , where T cut is set to 0.05. The fit of the 
numerical RGM data reproduces the analytical results of 
Ref. H3 for J 2 = 0, i.e., 6( J 2 = 0) = 2/?(J 2 = 0) = tt/2, 
with a precision of four digits. 

After having tested our fitting procedure by compar- 
ison with the analytical predictions for J 2 = 0, we now 
consider the frustrated model, where (to the best of our 
knowledge) no other results are available. From our nu- 
merical data for x an d £ we determine the J 2 dependence 
of the coefficients b and j3. We find that the numerical 
data for b and (3 obtained by the fitting procedure de- 
scribed above are very well described by a linear decrease 
in both parameters with increasing frustration, 

b = 2/3=-^(J 1 +2J 2 ). (11) 

Obviously, both parameters would be zero at the classi- 
cal transition point J 2 = 0.5, but they are still finite at 
the transition point Jf sa 0.44 of the quantum model. 
Hence, the exponential divergence is present in the full 
parameter range J 2 < Jf where the ground state is fer- 
romagnetic. Wc emphasize that this result is contrary 
to the behavior observed for the ID frustrated spin-1/2 
ferromagnet, where the critical properties change at the 
zero-temperature transition points We mention further 
that the leading coefficients ao and ao of the preexponen- 
tial functions, see the expressions for x and £ given above, 
vanish at the transition point Jf , whereas the next coeffi- 
cients a\ and ai remain finite at Jf . Hence, the preexpo- 
nential temperature dependence is changed approaching 
the zero-temperature transition. 

Let us mention, that the linear decrease in the coef- 
ficients b and j3, see Eq. pip , found by fitting the low- 
temperature behavior of x an d £ is the same as that ob- 
tained analytically for the zero-temperature spin-stiffness 
p s , see Sec. [TTJ This relation between p s and the diver- 
gence of the correlation length and the susceptibility is in 
accordance with general argument a 48 ' 49 concerning the 
low-temperature physics of low-dimensional Heisenberg 
ferromagnets. 

Another interesting quantity is the specific heat CV 
shown in Fig. [6] For J 2 = the specific heat ex- 
hibits a typical broad maximum at about T = 0.562. 
Increasing the frustration the height of this maximum 
becomes smaller and it is shifted to lower temperatures. 
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Interestingly, within our RGM approach the shape of the 
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FIG. 6: Specific heat calculated by RGM for J2 = 
0,0.1,0.2,0.3,0.35, and 0.4, from top to bottom. The inset 
shows the specific heat for J2 = 0.3 to 0.35 in steps of 0.01. 

Cv(T) curve changes for large frustration. For J2 > 0.34 
the specific heat shows two maxima, one at low tempera- 
tures (T < 0.193) and another one at high temperatures 
(T > 0.6). This extra low-temperature maximum signals 
the emergence of an additional low-energy scale when ap- 
proaching the transition point J|. Then many low- lying 
multiplcts appear above the fully polarized ferromagnetic 
ground-state multiplet. Hence, the appearance of the 
additional low-temperature maximum in Cy(T) can be 
attributed to a subtle interplay between all of these low- 
lying states. Note that a similar behavior was found for 
the ID caseM 



IV. SUMMARY 

In this paper we investigated the influence of the frus- 
trating NNN coupling J2 < Jf on the thermodynamic 
properties of the 2D spin-1/2 Heisenberg ferromagnet us- 
ing the spin-rotation-invariant Green's function method 
(RGM). The RGM estimate for the critical frustration 



J|, where the ferromagnetic ground state breaks down, 
is J 2 C w 0.44|Ji|. 

We tested the method by comparing RGM data for 
the spin-spin correlation functions and the uniform sus- 
ceptibility x calculated for finite lattices of N = 16 and 
N = 20 sites with corresponding exact-diagonalization 
data for the same lattice sizes and found a good agree- 
ment between both methods. However, the comparison 
of RGM data for finite lattices and for N — > 00 indicates 
strong finite-size effects at low temperatures. This leads 
to the conclusion that the finite systems of N = 16 and 
N = 20 are not representative for the low-temperature 
thermodynamics of large systems. 

As it is known from the Mcrmin- Wagner theorem^ 
the thermal fluctuations are strong enough to suppress 
magnetic LRO for the Heisenberg ferromagnet in dimen- 
sion D < 3 at any finite temperature. Due to frustration 
the fluctuations are further enhanced. Thus, frustration 
leads to a significant suppression of magnetic correlations 
at finite temperatures even if the ground state remains 
ferromagnetic. 

The low-temperature behavior of the susceptibility x 
and the correlation length £ in the ferromagnetic ground- 
state region J2 < Jf exhibits the exponential divergences 

x oc exp( " (|Jl j- 2j2) ) and £ cx cxp ( " (|Jl j~ 2J2) ). Al- 

though the numerator in the exponent decreases with 
growing J2 the exponential divergence perpetuates in the 
whole parameter region J 2 < J%- 

The specific heat calculated within the RGM exhibits 
a double-maximum structure for J 2 > 0.34 analogous to 
the ID model. 

To test our theoretical predictions it would be inter- 
esting to reconsider experimentally quasi-2D spin-1/2 
ferromagnets, such as K2CUF4, CS2CUF4, Cs2AgF4, 
La 2 BaCu05, and Rb2CrCLp22r— with respect to a 
possible frustration and its influence on low-temperature 
thermodynamics. 
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